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Abstract. We present time-series spectroscopy and multi-band photometry of the sdBV PG 1605+072 carried 
out simultaneously at the Calar Alto 2.2m and 3.5m telescopes. The periodogram analysis of the radial velocity 
curves reveals three frequencies at 2.078, 2.756, and 1.985 mHz for }lp and at 2.076, 2.753, and 1.978 mHz 
for H 7 . The corresponding radial velocity amplitudes are 12.7, 8.0, and 7.9 km/s for H/3 and 14.3, 6.5, and 
7.2 km/s for H 7 . Furthermore, we found five frequencies that are present in all wavelength bands of the BUSCA 
photometer. The frequencies detected in the radial velocity curves are recovered by the photometric measurements. 
Moreover, additional frequencies were present in the periodograms which could not be identified in all four bands 
simultaneously. The comparison of the amplitudes presented here with previous results from radial velocity and 
photometric observations of PG 1605+072 shows a significant change or even switch in the power of the modes 
within short time scales, i. e. about one year. No changes in frequency were registered and the phases of the modes 
show no wavelength dependency within our multi-band photometry. 
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1. Introduction 

Subluminous B (sdB) stars dominate the populations of 
faint blue stars of our own Galaxy and are found in the 
disk (field sdBs) as well as in globular clusters (Moehler 
et al. I1997[) . Furthermore, these stars play an important 
role trying to explain the "UV upturn phenomenon" ob- 
served in elliptical galaxies and galaxy bulges (Greggio & 
Renzini 19901 HU^jl . According to observations with the 
Ultraviolet Imaging Telescope (Brown et al. 1997) and the 
Hubble Space Telescope (Brown et al. 2000) sdB stars are 
sufficiently numerous to be responsible for the excessive 
UV flux. 

It is generally accepted that sdB stars can be identi- 
fied with models of the extended Horizontal Branch (EHB) 
burning He in their core fHeber lH)8 6 SafFer et al. Q994). 
The hydrogen envelope surrounding the core of about half 
a solar mass is very thin (< 2 % by mass) and therefore in- 
ert. These EHB stars will continue their evolution directly 
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towards the white dwarf graveyard avoiding the AGB and 
planetary nebula phases (Dorman et al. [1993). 

How sdB stars evolve towards the EHB with effective 
temperatures of up to 40 000 K remains a puzzle. The star 
must have lost all but a tiny fraction of the hydrogen en- 
velope at the same time as the He core has attained the 
minimum mass (w0.5Mq) required for the He flash. This 
challenges every mass loss mechanism in aspects of timing 
and effectivity. Recent findings (Maxted et al. 2001 . SafFer 
et al. 120011 Heber et al. 2002) emphasize the significance 
of close binary evolution. 

Some of the sdB stars were recently found to exhibit 
rapid multi-periodic light variations (P « 80-600 s) of 
low amplitudes (a few mmag). They form a new class of 
pulsating stars named after the prototype EC 14026 stars 1 
(Kilkenny et al. 1997). Since then a relatively large number 
of new sdB pulsators has been discovered. 31 are known 
today fCharpinet 120011 Piccioni et al EZUoTJl Silvotti et al. 
2000|. The observed brightness variations are caused by 
radial and non-radial, low degree and low order acoustic 
pulsation modes. The pulsations in these stars are driven 
by an opacity bump due to Fe and other metallic species 
(Charpinet et al. I1997|) at a temperature of « 2 x 10 5 K 
in the sdB envelope. 

1 EC 14026 stars are now officially V361 Hya stars 
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Stellar pulsations allow a direct insight into the struc- 
ture of such stars and therefore into the evolutionary his- 
tory. The frequencies or periods of the pulsation modes 
probe the chemical stratification and the mass which oth- 
erwise are difficult or even impossible to determine. The 
power of asteroseismological tools has been demonstrated 
in the field of pulsating white dwarfs for which stellar pa- 
rameters like mass, luminosity or thickness of the enve- 
lope were derived (e. g. Winget et al. I1991[l . In the case 
of variable sdB stars these parameters will constrain the 
evolutionary history and consequently shed more light on 
the origin of these stars. 

Identification of pulsation modes (characterized by 
spherical harmonics with the indices / and m) is a pre- 
requisite for asteroseismology. Brassard et el. (|2001[l have 
successfully carried out an asteroseismological analysis for 
PG 0014+067. For the first time they were able to deter- 
mine the stellar mass (M*/M = 0.490 + 0.019) as well as 
the envelope mass (log (M onv /M ) = -4.31 ± 0.22) and 
both are in excellent agreement with predictions from evo- 
lutionary models (Dorman et al. 1993). 

Pulsations produce not only photometric variations 
but also line profile variations that offer an alternative ap- 
proach towards mode identification. PG 1605+072 is the 
ideal target for this application: it has the longest pul- 
sation periods (« 500 s) which enables spectra with rea- 
sonable S /N to be obtained within each pulsation period. 
Moreover, this star has the largest variations of all known 
sdBVs (0.2 mag in the optical) and by far the richest fre- 
quency spectrum (> 50 modes, Kilkenny et al. I1999f) . A 
recent spectroscopic study revealed this star to be a rather 
rapid rotator (vsini = 39 km/s, Heber et al. H999|l which 
will complicate the identification of modes due to non- 
linear effects on mode splitting. 

This work presented here serves as a feasibility study 
in order to find out whether an asteroseismologic anal- 
ysis of PG 1605+072 is possible. For this purpose we 
have done simultaneously time-series spectroscopy and 
multi-band photometry (observations and reductions in 
Sect. |2l . The analysis of the results is presented in Sect.|3] 
Previously, other groups have done radial velocity studies 
(O'Toole et al. lMlfHIM^l Woolf et al. lMHSfl or photomet- 
ric campaigns (e. g. Kilkenny et al. H999|) on PG 1605+072. 
Simultaneous multi-band photometry has not been ob- 
served before. This enables us to study the temporal evo- 
lution of the frequencies and amplitudes of the pulsation 
modes (Sect.0J. Finally, we discuss our results and give a 
brief outlook to future work in Sect. 

2. Observations and data reduction 

2.1. Time-series spectroscopy 

We obtained time resolved longslit spectra for the sdBV 
PG 1605+072 during one observing run. The observations 
were carried out with the 3.5m telescope at the DSAZ, 
Calar Alto, Spain and the TWIN spectrograph (two SITe- 
CCDs with 2048 x 800 pixels a 15 /im). The time alloca- 



tion committee awarded two nights on May 14 and 15, 
2001. Because of rain during the first night data were only 
obtained during the second night (5h27m of data). The 
passage of clouds and variable seeing affected the quality 
of the data considerably. In order to achieve an adequate 
time resolution observations were carried out in trailing 
mode which means that the telescope tracking is tuned 
such that the star moves slowly along the slit in N-S di- 
rection. The drift velocity was chosen at 270" per hour 
so that a time resolution of 15 s per pixel resulted. The 
time resolution is dependent on the seeing conditions and 
thus varied during the night. The slit width was 1.5". We 
chose the gratings T05 and T06 (36 A/ mm) correspond- 
ing to a spectral resolution of « 1 A. The spectral ranges 
from 3860 - 4960 A in the blue and from 5880 - 6920 A in 
the red were covered. Moreover, to decrease the readout 
time and the noise level the data were binned (1x2). 

2.2. Multi-band photometry 

Simultaneously with the time resolved spectroscopy of the 
target star its light curve was measured in four wave- 
length bands at the 2.2m telescope at Calar Alto. At this 
telescope five nights between May 14 and 18, 2001 were 
awarded. No data could be taken during the first night and 
from the middle of the fourth night until the end of the run 
due to bad weather conditions. 12 h 27 m of pure observa- 
tion time resulted in 880 measurements. The telescope was 
equipped with the new multi-band photometer BUSCA 
(Reif et al. 1999, Cordes et al. in prep.) which is able to 
measure in four photometric bands simultaneously. Three 
beamsplitters split the incoming light into four beams that 
feed four 4kx4k CCDs. Only Str0mgren filters were avail- 
able for BUSCA. Use of these filters would have further re- 
duced the S/N. Therefore, no filters were inserted and the 
beamsplitters served as non-overlapping broad band fil- 
ters. The transmission curves of the four bands are shown 
in Fig. These curves convoluted with the quantum effi- 
ciency curves of the CCDs provide the pass bands for our 
observations. The four wavelength bands are denoted by 
W B \ 'B B ', 'i? B ' and 'NIR B ' throughout the paper fol- 
lowing the notation of Cordes et al. (in prep.). We chose 
an exposure time of 15 s which provides a good S/N and a 
reasonable coverage of the main frequencies. 2x2 binning 
(resolution: 0.17" per unbinncd pixel) and the windowing 
option (reading out 346x400 pixels of the CCDs) helped 
to reduce the full cycle time to 51 s. 

The brightness variations of PG 1605+072 were de- 
tected by relative photometry. For this purpose two stars 
of similar brightnesses in the R band (Ri = 13.2 and 
R2 = 13.4) next to our target star were used. 

2.3. Data reduction 

Data reduction was done using the IDL based software 
SPEX (Long slit SPectrum Extraction package 2 ), and 

2 see http://astro.uni-tuebingen.de/~schuh/spex/index.html 
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Fig. 1. Transmission curves of the four BUSCA wave- 
length bands. The four ranges are denoted by 'UVg' (cen- 
tral wavelength: 3 600 A), 'B B ' (4800 A), 'i? B ' (6 300 A) 
and WJi?g' (8 000 A) although these are not comparable 
to the classical Johnson UBVRI system. 

TRIPP (Time Resolved Imaging Photometry Package, 
see Schuh et al. lT535|l . TRIPP is based on the C CD pho- 
tometry routines written by R.D. Geckeler l|1998l) and per- 
forms aperture photometry. The most important step is 
the determination of the relative flux of the target star 
with one or more stars as comparison objects. The de- 
tection of variations of the order of a few mmag is only 
possible when comparison stars and sky background are 
recorded simultaneously (especially under comparatively 
poor conditions). 

SPEX is designed for the reduction of long slit spec- 
tra. The reduction comprises of the standard procedures 
flatfielding, biasing, cosmic ray event extraction, and 2-D 
wavelength calibration. The latter is very important be- 
cause we extract spectra from single CCD rows in order 
to establish the time-series. Therefore, we have to ensure 
that the position of the comparison lines perpendicular 
to the dispersion direction doesn't change. In fact, these 
positions vary and thus we correct for that by fitting a 
polynomial. After the calibration of one CCD row the re- 
sulting dispersion relation is applied to every single CCD 
row. 

One crucial point was the calibration of the time axis. 
As a consequence of the trailing mode observation the 
seeing disk moves slowly across the CCD and therefore 



the signal smears over the trailing direction. The initial 
point of time of the measurement and the cycle time were 
determined by fitting Gaussian profiles to the rising flank 
of the star's signal. Finally, radial velocities were derived 
by fitting Lorentz profiles to the Balmer lines H Q -Hs as 
well as to He i 4471 A and He n 4686 A. 

3. Analysis of the time-series 

Both the spectroscopic as well as the photometric time- 
series were analysed with the program package TRIPP. It 
enables the calculation of periodograms, confidence levels 
and fits with multiple sine functions (see Dreizler et al. 
2002 for more details). 

3.1. Spectroscopy 

We were able to derive accurate radial velocity curves for 
the Balmer lines H^ and H 7 (see Fig. EJ. The frequency 
resolution for this run is calculated to be 51 /iHz. Trailing 
mode observations require stable weather conditions so 
that the intensity of the star's signal is almost constant 
during one exposure. Thus bad seeing and transparency 
changes due to the passage of small clouds influence the 
quality of the data rather strongly. Because of these dis- 
turbances the quality of the radial velocity curves varies 
with time. The radial velocity curves of the other Balmer 
lines and the He lines turned out to be too noisy for a 
quantitative analysis and are not further discussed. Better 
conditions should allow a more extensive study of radial 
velocity changes in this star. 

The Lomb-Scargle periodograms calculated from the 
radial velocity curves (see Fig. show peaks clearly sur- 
passing the 3er confidence level. Applying a prewhitening 
procedure to the data three frequencies and their ampli- 
tudes (see Tabled were found for both Balmer lines in- 
vestigated. Other peaks that seem to lie above the 3<r level 
could not be extracted by our procedure. The dominant 
frequency at 2.076 mHz has the largest velocity ampli- 
tude with 12.7 km/s for Hg and 14.3 km/s for H 7 . The 
amplitudes of the other two frequencies are lower than 
the strongest one with values ranging from 6 to 8 km/s. 
We determined the amplitude accuracy AA by calculating 
the median value of the white noise in the frequency range 
3-7 mHz where almost no power arises in the frequency 
spectrum. According to this, we found AA= 1.5 km/s for 
H/3 and AA = 1.0 km/s for H 7 . 

3.2. Multi-band photometry 

BUSCA is a unique instrument which enables the mea- 
surement in four different wave bands simultaneously. As 
a result we obtained four light curves. 

Fig. 0| shows the Lomb-Scargle periodograms of all 
BUSCA bands. They are quite similar to the periodograms 
derived from the radial velocity curves (Fig. |3J. This 
data set spans over three nights so that daily aliases are 
clearly visible and the frequency resolution is much better 
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Fig. 2. Radial velocity curve for H^ (upper panel) and H 7 
(lower panel); artefacts due to poor weather conditions 
were not removed 

at Ai> — 5.68 /iHz (compared to 51 /zHz for the spec- 
troscopy). Again we applied the prewhitening technique 
in order to remove all significant peaks from the peri- 
odogram and to obtain the amplitudes for each frequency. 
As before, the horizontal line in the diagrams represents 
the 3cr confidence level above which we assume the de- 
tected frequencies to be real. In all four wave bands five 
peaks with the same frequencies can be identified (see 
Tabic [5] and The dominant frequency is again found at 
2.076 mHz and therefore confirms the results from spec- 
troscopy. Furthermore, additional frequencies were found 
in the region around 2.74-2.78 mHz but these peaks are 
closely spaced so that a corresponding identification in all 
BUSCA bands due to the medium frequency resolution 
was not possible. Peaks that fall below the 3cr confidence 
level were not removed from the periodograms. 

The amplitudes of the brightness variations are 
measured in fractional intensity. The light curves are 
normalized to the fraction of intensity <5I/I and the 
amplitudes were then converted to mmag. This 
unit will be used throughout the whole paper. The 
accuracy of these amplitudes are calculated in the same 
way as we did it for the radial velocity amplitudes. Here 
we used the median level of the white noise in the range 
3-5 mHz. The accuracies for the BUSCA wavebands are 
1.52mmag for 'C/Vg', 1.53 mmag for 'B-g\ 1.12mmag for 
'i?g' and 1.37 mmag for , NIRq\ respectively. 



Fig. 3. Lomb-Scargle periodogram of the radial velocity 
curve of H^ (upper panel) and H 7 (lower panel). The 
power spectral density (psd) is a measure for the probabil- 
ity that a period is present in the radial velocity curve or 
light curve. The horizontal line represents the confidence 
level of 99 % (3<r-level). 



Fig. shows the semi amplitudes of four selected fre- 
quencies as a function of effective wavelength of the bands. 
In Fig. El we display the relative change of the semi am- 
plitude of each waveband. The deviation with respect to 
the mean is largest for the , UVq band. The other chan- 
nels, considered separately, behave rather similar showing 
much smaller deviations from the mean brightness. This 
is explained through the fact that the 'UVq band lies 
blueward to the Balmer jump and the other redward of 
it. The opacity changes a lot across this wavelength range 
and thus the stellar flux originates from different atmo- 
spheric depths. 

Furthermore, we used the phases which are delivered 
by the sine fit procedure (see Table[2]and[3J) in order to test 
whether there is any wavelength dependency. The phase 
values are normalized to unity. Fig.0shows the deviations 
of the phases with respect to the mean value of all four 
wavebands for the four selected frequencies of Fig. [3] The 
error bars are calculated from 

A = Re{A} sin0. (1) 



S. Falter et al.: Simultaneous time-series spectroscopy and multi-band photometry of the sdBV PG 1605+072 



5 



The phase error A(f> can be calculated from the amplitude 
uncertainty AA according to 

A0 = arctan ■ ( 2 ) 

The two strongest frequencies at 2.076 and 1.986 mHz 
show no wavelength dependent behavior within the er- 
ror margins. The deviations for the other two frequencies 
from the mean value are slightly larger than the formal 
errors. However, we regard this as insignificant, and con- 
clude that we cannot find any wavelength dependency of 
the phases of the modes. 

4. Comparison with previous investigations 

Other groups have carried out similar studies already, ei- 
ther time-series spectroscopy (O'Toole et al. 120001 O'Toole 
et al.EOOS Woolf et al. or photometry (Kilkenny et 

al.EHnU- O'Toole et al. fZUHty have collected monochro- 
matic photometric data, too, but they are not strictly si- 
multaneous. We shall compare our results with their re- 
sults in this section. 

The photometric multi-site campaign of Kilkenny et 
al. I|1999f) discovered more than 50 frequencies in their 
frequency analysis of PG 1605+072. The five frequencies 
(2.076, 1.9865, 2.102, 2.2695 and 1.8912mHz) detected in 
all BUSCA bands of our data were already present in 
their data. Some of the others which we found addition- 
ally (e.g. 4.0631 mHz) can be identified with frequencies 
from Kilkenny et al. Ijl999(l and others can not (in par- 
ticular at 2.7866 mHz). A closer inspection of these 
frequencies shows that all of them are very close 
to one-day aliases (±11. 57 /iHz) of the extracted 
frequencies of our or of Kilkenny's data to within 
the frequency resolution of 5.68 /iHz. Small devia- 
tions of our measurements from those of Kilkenny et al. 
(1999) can be explained by the shorter time spanned by 
our data (« 48 h compared to « 217.6 h, frequency resolu- 
tion: 5.68 /iHz compared to 1.28 /iHz). In contrast to our 
analysis, Kilkenny et al. I|1999|) were able to resolve closely 
spaced frequencies, e.g. at 2.1017/iHz and 2.1033 /iHz (see 
Table Since 1997 the relative power within one blend 
of frequencies may have changed so that the frequency of 
unresolved peaks in the periodogram changes slightly be- 
cause the contribution of the weaker feature at that time 
may have grown in the meantime. We conclude that the 
frequencies measured in the BUSCA light curves are con- 
sistent with those of Kilkenny et al. <|1999p . The frequen- 
cies are stable in time to within our measurement errors. 

Also the frequencies discovered in the radial veloc- 
ity curves are consistent with both the BUSCA and the 
Kilkenny values. The time basis of the spectroscopic data 
is too short in comparison to the long-term photometry 
to resolve the frequencies. This is also true for most of 
the other spectroscopic feasibility studies mentioned be- 
fore. Only O'Toole et al. (2002) had a sufficiently long 
baseline of « 70 d. Thus, we discuss only the relative dis- 
tribution of the power of the pulsation frequencies. As 



can be seen in Fig. [3] the main power is located around 
2.076 mHz. Significant power arises in the frequency range 
2.74-2.78mHz but it falls off compared to the former. 
The same distribution of pulsation power was detected 
by the multi-site photometric campaign of Kilkenny et al. 
( 199SJ) who carried out their observations in April and 
May 1997. Two years later, in July and August 1999, 
O'Toole et al. <|2000f) recovered the same frequency pattern 
but discovered that the amplitudes had changed and, in 
particular, that no power was detected at 2.076 mHz. 
Another year later (in May 2000) O'Toole et al. l|20Tj2jl 
and Woolf, Jeffery & Pollacco (|2)02) gathered data again. 
Both groups discovered a drastic change in the relative 
power distribution: The power at 2.076 mHz appeared 
again and was at that time weaker than that in the fre- 
quency range 2.74 -2.78 mHz. This means that the power 
switches within a few years. Our data which were taken 
in May 2001, i. e. another year later, show that the power 
distribution switched back to that of 1997 measured by 
Kilkenny et al. l|F9"9"9"jl . 

5. Conclusion 

Our spectroscopic feasibility study based on data taken 
in May 2001 confirms earlier studies from observations in 
July/August 1999 and May 2000. All these investigations 
have shown that it is possible to measure accurate radial 
velocity curves for PG 1605+072. Three to five frequen- 
cies could be detected in the periodograms. The detected 
frequencies found in the radial velocity studies are very 
stable (see Table |21 and 0) . However, the amplitudes of 
the frequencies show drastic temporal changes over a few 
years. We also have simultaneously measured photomet- 
ric variations of PG 1605+072 in four wavebands for the 
first time. The amplitudes in the 'UV-q band are found to 
be significantly larger than at larger wavelengths ('-Bg', 
'i?g' and 'NIR-q'). No phase shifts between the four light 
curves can be detected. 

The time basis of all observations is much too short 
to resolve the full frequency spectrum of this star. More 
than 50 frequencies have been resolved in the long pho- 
tometric study of Kilkenny et al. (19921). Consequently, 
a new data set with a longer time basis is needed. For 
that reason we organised a combined photometric and 
spectroscopic multi-site campaign (PI: U. Heber), termed 
the MSST (Multi Site Spectroscopic Telescope) project 
(Heber et al. 2003 i. These observations covered a period 
of one month in May/ June 2002. This project is unique 
in that it combines photometric and spectroscopic data. 
The data still have to be reduced and analysed. A lot of 
work remains to be done. Hopefully, we will be able to 
identify pulsation modes within this large data set and 
carry out a detailed asteroseismological investigation of 
PG 1605+072. The feasibility study presented in this pa- 
per in concert with all other equivalent studies mentioned 
in Section 01 have paved the way for this aim. 

Acknowledgements. The authors would like to thank the 
Calar Alto staff for their support during our obser- 
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Table 1. Radial velocities (RV) of PG 1605+072 derived for and H 7 and for comparison of the other radial velocity 
studies. Woolf, Jeffery & Pollacco l|2002p measured their radial velocities from shifts of the whole spectrum. The radial 
velocities of O'Toole et al. I|2002|) are measured for all Balmer lines (velocity error: w 0.4km/s). 
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Table 2. Frequencies and corresponding periods as well as amplitudes for PG 1605+072 in the first two BUSCA 
wavebands ('C/Vg' and '-Bp,')- The errors of the periods and amplitudes are the formal fit errors from the sine fit 
procedure. The left column shows the corresponding values derived in Kilkenny et al. I|1999f) as a comparison to the 
values found in this work. The two frequencies marked with a star can not be resolved in our data but can in the data 
of Kilkenny (see Sect. 0] for the discussion). The phases derived within the sine fit procedure are denoted by 4> and are 
discussed in Sect. [3] as well as the determination of the phase errors. 
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Table 3. Frequencies and corresponding periods as well as amplitudes for PG 1605+072 in the last two BUSCA 
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Fig. 5. Semi amplitudes of four frequencies as a function 
of effective wavelength. The error bars are la errors and 
are calculated by means of a x 2 method. 
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Fig. 4. Lomb-Scargle periodogram of the light curves of 
the 'J/Vb', 'B b \ 'R b ' and W/i? B ' band (from top to 
bottom) of the BUSCA camera. The horizontal line in 
each panel represents the confidence level of 99 % (3a- 
level) . 
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Fig. 6. Relative semi amplitudes of the four frequencies 
from Fig. [5] These are calculated by Ay = (y — y) /y. y is 
the semi amplitude of the wavebands. The error bars are 
the formal la errors of the sine fit procedure. 



Fig. 7. Deviation of the phases from the mean value cal- 
culated from the four frequency bands for the frequencies 
given in Fig. The error bars are calculated from the 
amplitude errors as described in the text. 
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